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We determine the quark distributions and structure functions for both unpolarized and polar- 
ized DIS of leptons on nucleons and nuclei. The scalar and vector mean fields in the nucleus modify 
LT* ■ the motion of the quarks inside the nucleons. By taking into account this medium modification, 

we are able to reproduce the experimental data on the unpolarized EMC effect, and to make pre- 
dictions for the polarized EMC effect. We discuss examples of nuclei where the polarized EMC 
^ ' effect could be measured. We finally present an extension of our model to describe fragmentation 

H . functions. 



1 Introduction 

In order to describe nonperturbative effects like spontaneous chiral symmetry breaking and nuclear 
binding on the level of quarks, effective chiral quark theories are powerful tools. A prominent example 
is the EMC effect, which has clearly shown that the quark distributions of bound nucleons differ from 
those of free nucleons [1] . It has been shown recently ^ that this effect can be explained if one takes 
into account the response of the quark wave function to the nuclear environment, that is, to the nuclear 
mean fields, and that the same mechanism gives rise also to medium modifications of the polarized quark 
distributions and structure functions. In this work, we will focus on the predictions for the polarized 
EMC effect, and briefly discuss extensions to describe transversity distributions |3j and fragmentation 
functions \^. 
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2 Quark distributions and structure functions 

In this work we will mainly be concerned with the following EMC ratios: 






Here x is the Bjorken variable for the nucleon, and xa is A times the Bjorken variable for the nucleus of 
mass number A, so that < xa < A. The unpolarized and polarized structure functions of the nucleon 
are denoted as F2a and gia respectively (a = p,n), while F2A and (^f^ ^^^ the corresponding structure 
functions of the nucleus A with spin projection H = —J, ■ ■ ■ + J along the direction of the incoming 
electron. The polarization factors of protons and neutrons, which appear in the denominator of the 
spin dependent EMC ratio, are defined as twice the expectation values of the proton and neutron spin 
operators between the polarized nuclear states. Both ratios in Eq. ([1]) are defined so that they become 
unity in a naive single particle model based on nonrelativistic nucleons. 

The parton model expressions for the nuclear structure functions are very similar to those of the 
nucleon fS], for example 

9ia{xa) = ij^ejAgf (xa) = ^^24 {q^^{xA) - q^^ixA)) • (2) 

<? g 

Here qA-ii^A) is the probability to find a quark with momentum fraction Xa/A and Sz = 1/2 in the 
nucleus A with J^ = H, and similar for q^^{xA)- 

It is important to keep in mind the following two points: First, usually only a few valence nucleons 
(or holes) contribute to the nuclear polarization, and therefore ^^f^ is of order 1/A relative to F2A, where 
all nucleons contribute. Second, the structure function of a free proton is larger and much better known 
than the neutron structure function. Therefore, possible candidates for the observation of the polarized 
EMC effect are stable nuclei which are not too heavy, and where the polarization is dominated mainly 
by the protons. 

Nuclear spin sums are interesting quantities, which have not yet been explored in detail. The 
isoscalar and isovector combinations are 

f dxA (AuHxa) + AdiixA)) = S (P/ + P„^) , 
J dxA {AuHxa) - AdiixA)) = QA (P/ - Pil) , (3) 

where we assumed for simplicity that only a single nucleon state contributes to the nuclear polarization. 
The first relation contains information on the quark spin sum in a bound nucleon (S), and the second one 
on the axial coupling constant of a bound nucleon {qa)- The latter is related to nuclear Gamow- Teller 
matrix elements [6], and establishes an important link between quark physics and nuclear structure 
physics. 

3 Model calculations 

We describe the nucleon as a bound state of a quark and a diquark by using the Faddeev framework in 
the Nambu-Jona-Lasinio (NJL) model. Although the relativistic Faddeev equation, which is represented 
graphically in Fig. 1, can be solved exactly in the NJL model [7], for our applications to nuclei we limit 
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Figure 1: Graphical representation of the Faddeev equation. The single line denotes a quark, 
and the bold line a diquark. 
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Figure 2: Feynman diagrams for the quark distribution functions in the nucleon. The 
operators inserted into the quark lines are explained in the text. 



ourselves to the static approximation, where the momentum dependence of the quark exchange kernel is 
neglected |8]. We take into account the scalar and axial- vector diquark channels, and avoid unphysical 
quark decay thresholds by introducing an infrared cut-off, in addition to the ultraviolet one, in the 
proper-time regularization scheme [9]. 

By using the quark-diquark vertex functions, we evaluate the Feynman diagrams shown in Fig. 2 
with the operator insertions (7_, 7.75, 7.7^75) x b{x — k^jpJ), to get the unpolarized, polarized, and 
transversity quark distributions, respectively. Here k and p are the momenta of the quark and the 
nucleon, and a_ denotes the light-cone minus-component of a 4-vector a. It is straight forward to 
extend this calculation to the case of a bound nucleon by including the mean nuclear scalar and vector 
fields into the quark propagators of Fig. 2. 

Finite nuclei in our present approach are described in a simple independent particle picture: We 
assume Woods-Saxon scalar and vector potentials for nucleons, with depth parameters determined 
from our earlier self-consistent nuclear matter calculations [n], and standard values for the range and 
diffuseness parameters [10]. After solving the Dirac equation with these potentials, we calculate the 
expectation values of the scalar and vector potentials for each nucleon orbit, and use the quark-diquark 
(Faddeev) equation to translate them into the average scalar and vector fields for quarks. These 
average fields are then used in the quark propagators of Fig. 2. Finally, we calculate also the light-cone 
momentum distributions of the nucleons [2] , and obtain the quark distributions in the nucleus by using 
the convolution formalism. 

The parameters of the model are determined as usual from the properties of the pion and the free 
nucleon. In particular, the 4-Fermi coupling constants in the scalar and axial-vector diquark channels 
are fitted to the mass and the axial- vector coupling constant of a free nucleon [TT]. Then there is 
only one free parameter in the calculation of the quark distribution functions, which is the model scale 
(Qo) needed to perform the Q^ evolution [J. By fixing this scale to the same value as our constituent 
quark mass (400 MeV), we obtain a very good description of the empirical quark distributions in the 
free nucleon. This is shown in Figs. 3 and 4, where the results after the Q^ evolution (solid lines) are 
compared to the parametrizations of Refs. [121 [H] (dashed lines). Recently obtained results for the 
transversity distributions are shown in Ref. [3]. The results are similar to the helicity distributions. 



^We use the computer code of Ref. [?] to perform the Q^ evolution in NLO. 
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Figure 3: Unpolarized valence quark distributions in the free proton. 
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Figure 4: Polarized valence u-distributions (upper part) and ^-distributions (lower part) in 
the free proton. 



which is contrary to the analysis of Ref. [15]. Further investigations on this point are necessary. 



4 Results for nuclear quark distributions and structure func- 
tions 

In this section we will show our results for the medium modifications of unpolarized and polarized quark 
distributions and structure functions. 

Fig. 5 shows the unpolarized valence up quark distribution in the nucleus ^^B, and Fig. 6 shows 
the polarized up and down quark distributions for the same nucleuso The dotted hues are the free 
distributions, i.e., the results obtained by neglecting Fermi motion and medium effects. The dash- 
dotted lines include the effect of the scalar mean field, the dashed lines include further the Fermi 
motion, and finally the solid lines incorporate also the effect of the vector mean field. 



^Thc distributions and structure functions shown in Figs. 5-9 refer to the leading niultipolcs {K = for the unpolarized, 
K = 1 for the polarized case), which are linear combinations of the corresponding quantities in the helicity (H) basis. 



For details, see Ref. [S]. 
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Figure 5: Unpolarized valence u-distribution in B 
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Figure 6: Polarized valence w-distribution (upper part) and d-distribution (lower part) in 
"B. 



Comparing the dotted (free result) and the solid (full result) lines in Fig. 5, we see that the unpo- 
larized distribution becomes softer in the nucleus, and that the vector potential plays a very important 
role to describe this shift to smaller x [TB1|T7]. The main features shown in Fig. 5, namely a quenching 
of the distribution function at large x and a small enhancement at smaller x, are consistent with the 
EMC effect. On the other hand, Fig. 6 shows that the polarized quark distributions are quenched in 
the nucleus for all values of x, which implies a reduction of the quark spin sum for a bound nucleon 
compared to the free nucleon case. In other words, in the medium a part of the quark spin is converted 
into orbital angular momentum. 

The resulting EMC ratios of Eq. ([1]) for ^^B are shown in Fig. 7. It is seen that the polarized EMC 
effect is predicted to be larger than the unpolarized one. As further possible candidates to measure the 
polarized EMC effect, we show the results for ^Li and ^^Al in Figs. 8 and 9. We see that the difference 
between the two EMC ratios becomes more pronounced as the mass number increases. 

The spin sums for the free and the bound nucleon are listed in Table 1. Note that the quantities in 
the second to fifth columns are defined by dividing out the nuclear polarization factors from the nuclear 
spin sums (see Eq. (|3])), and therefore directly reflect the medium modifications. The last row shows 
the limit of infinite nuclear matter. We see that the quark spin sums are appreciably quenched in the 
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Figure 7: EMC ratios for B. The experimental data refer to C 
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Figure 8: EMC ratios for ^^Al 



medium. The last two columns of Table 1 show the tensor charges for the free nucleon [3] and a nucleon 
bound in infinite nuclear matter. 



5 Extension to fragmentation functions 

Here we wish to discuss a framework to extend the model to fragmentation functions. Numerical results 
for fragmentation functions will be presented in a future publication [1]. 

There exists a relation between the quark distribution inside a hadron {fq{x)) and the quark frag- 
mentation function into a hadron {D^[z)), which is known as the Drell-Levy-Yan relation y^. Let us 
discuss this relation, starting directly from the operator definitions 



fii^) = 2 5Z„ ^ (P-^ -P-+ Pn-) (pl-^lPn) 7^ (Pnl^b) , (4) 

Here \p) denotes the hadron state (we assume a nucleon for definiteness) , and ^^ is the sum over 
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Figure 9: EMC ratios for ^Li. The experimental data refer to ^Be. 



Table 1: Quark spin sums for a free proton and a proton bound in the nuclear medium. 





Am 


Ad 


S 


9A 


Atu 


Ard 


p 


0.97 


-0.30 


0.67 


1.27 


1.04 


-0.24 


^Li 


0.91 


-0.29 


0.62 


1.19 






iiR 


0.88 


-0.28 


0.60 


1.16 






15N 


0.87 


-0.28 


0.59 


1.15 






2^1 


0.87 


-0.28 


0.59 


1.15 






nucl. matt. 


0.74 


-0.25 


0.49 


0.99 


0.93 


-0.23 



the intermediate states |p„), including an integral over the momentum and sums over spin and isospin 
projections. For later convenience, the sum in Eq. ([5]) is taken over the antiparticle states (p^). 
We can express the matrix element in the distribution function (jlj) as 



(Pnl'iPlp) = ^iP,Pn) ^/%Un{Ps) 



(6) 



where the Dirac matrix r(p, p„) is the Fourier transform of the Green function (p„ | (^(0)<l>(2;)) | O) 
with $ the nucleon field, and y^A^ is a normalization factor for the nucleon spinor u^. Using crossing 
and charge conjugation symmetries, the matrix element in the fragmentation function ^ can then be 
expressed as 



(oiV'b,^;) 



NpVN{ps)T{-p,Pn) C 



(7) 



where C = i'j^'jo. Inserting these matrix elements and their conjugates into (jl]) and (I5l), it is easy to 
verify that 
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where p -^ —p means to reverse all 4 components of p^, and after this replacement p^ = E^ip) > 0. 
The effect of the replacement p -^ —p on the distribution function Eq. (jl]) is seen most easily by 
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expressing it in terms of F and the quark momentum k = p — pn as follows: 

X 5(A;_ -p_x) Tr {(^ + Mn)T{p,p - k)-f+T{p,p - k)) , 
= e(l-a;)F(x). (9) 

Here we expressed the integrand using light-cone variables, and the "energies" of the intermediate 
state (invariant mass M„) and the nucleon are defined by e„(p„) = (p^_|_ + M^)/(2p„_) and eN{p) = 
(pi + Mjf)/{2p^). From ([8]) and iQ we obtain finally 

D^^{z) = -e{l-z)'-F(^x = ^^ . (10) 

(For spin zero bosons, there is no minus sign in Eq. (ITOl).) This shows that fj} and D^ are essentially 
one and the same function, defined in different regions of the variable. This important result allows us 
to extend our investigations on the distribution functions presented in this paper to the fragmentation 
functions. The numerical results and detailed discussions will be presented in a future publication [4J. 
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